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We report the result of a search for neutrino oscillations using precise measurements of the recoil 
electron energy spectrum and zenith angle variations of the solar neutrino flux from 1258 days of 
neutrino-electron scattering data in Super-Kamiokande. The absence of significant zenith angle 
variation and spectrum distortion places strong constraints on neutrino mixing and mass difference 
in a flux-independent way. Using the Super-Kamiokande flux measurement in addition, two allowed 
regions at large mixing are found. 



For over 30 years, measurements of the solar neutrino 
flux [0, ^, H, ^, 1^ have been significantly below the pre- 
diction of the Standard Solar Models (SSMs)§, §. Neu- 
trino flavor oscillations, similar to those seen in atmo- 
spheric neutrinos 0], are a natural explanation for this 



discrepancy. This type of flavor conversion is inher- 
ently energy-dependent. Since Super-Kamiokande (SK) 
measures the energy of the recoil electron from elastic 
electron-neutrino scattering, it has sensitivity to this en- 
ergy dependence. In addition to a conversion in vacuum. 
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TABLE L Flux, uncertainty and definition of zenith angle and energy bins. The systematic uncertainty in the last two columns 
is split into energy-uncorrelated and energy-correlated uncertainty. The systematic uncertainty is assumed to be fully correlated 
in zenith angle. 



Fluxzbstatistical uncertainty in units of SSM 



syst. uncert. in % 



cos^z-Range 



Day Mantle 1 Mantle 2 Mantle 3 

-1.00-0.00 0.00-0.16 0.16-0.33 0.33-0.50 



Mantle 4 Mantle 5 Core 
0.50-0.67 0.67-0.84 0.84-1.00 



Energy- 
uncorr. correlated 



5.0-5.5 MeV 
5.5-6.5 MeV 
6.5-8.0 MeV 
8.0-9.5 MeV 
9.5-11.5 MeV 
11.5-13.5 MeV 
13.5-16.0 MeV 
16.0-20.0 MeV 



0.436±0.046 

0.431±0.022 0.464±0.060 0.410±0.055 0.442±0.048 
0.461±0.013 0.524±0.036 0.506±0.033 0.438±0.028 
0.437±0.014 0.449±0.038 0.482±0.036 0.460±0.031 
0.434±0.015 0.432±0.042 0.493±0.040 0.446±0.034 
0.456±0.026 0.496±0.071 0.290±0.055 0.394±0.053 
0.482±0.056 0.532±0.155 0.775±0.171 0.685±0.141 

0.476±0.149 



0.453±0.048 0.495±0.054 0.434±0.058 
0.466±0.027 0.424±0.030 0.409±0.033 
0.503±0.031 0.461±0.034 0.439±0.037 
0.448±0.034 0.435±0.037 0.484±0.044 

0.477±0.056 0.439±0.061 0.465±0.068 
0.607±0.130 0.471±0.128 0.539±0.153 



+3.9 
-3.1 

-fl.5 
-1.4 

±1.4 

±1.4 

±1.4 

±1.4 

±1.4 

±1.4 



+0.25 
-0.21 
+0.30 
-0.26 
+0.77 
-0.75 

±1.6 

+3.1 
-2.9 
+5.5 
-5.0 
+9.2 
-8.3 
+16 
-14 



a matter-induced resonance in the sun||9| may sufficiently 
enhance the disappearance probability of solar neutri- 
nos even for small neutrino mixing. For some oscilla- 
tion parameters, matter-enhanced oscillations within the 
earth can lead to a different flux during day-time than 
during the night. Therefore, either a distorted energy- 
dependence or a zenith angle variation of the solar neu- 
trino flux would be considered evidence for oscillation. 
In addition, the shape of the distortion and zenith an- 
gle variation would determine the oscillation parameters 
independently from uncertainties in the SSM's flux pre- 
diction. As a real-time solar neutrino experiment, SK 
can study the solar zenith angle flux variation. 

Super-Kamiokande started taking data in April, 1996. 
It has since confirmed the deficit of solar neutrinos 
measured the recoil energy spectrum and carried out 
an initial search for zenith angle variation ||ll[| . In this 
report, we analyze the first 1258 days of data (May 31^^, 

1996 through October 6^^, 2000) using both spectral dis- 
tortion and zenith angle variation. The total number of 
solar neutrino events above a threshold of 5 MeV of recoil 
electron energy is IS'idi'^f^Q. The resulting flux of ®B so- 
lar neutrinos, (2.32j:° °^) x 10V(cm2s)[|, is OAbl^l^H 
of the flux predicted by the reference SSM(BP2000)@1. 

The sample is divided into seven zenith angle bins (one 
day bin and six bins in cos9z for the night); within each 
zenith angle bin, the data are divided into eight recoil 
electron bins. We will refer to this binning of the data 
as the "zenith angle spectrum" (see Fig. p. We define 
the zenith angle 9z of an event as the angle between the 
vertical direction and the solar direction at the time of 
the event. Day events have cos6'^ < and night events 
cos^z > 0. The size of the sample (already divided into 
seven zenith angles) does not allow a subdivision into 19 
energy bins shown in [||. Due to this statistical limita- 
tion the lowest (5.0-5.5 MeV) and the highest (16.0-20.0 



MeV) energy bin combine the fiux of all zenith angles. 
Table | shows the flux, statistical and systematic uncer- 
tainty for all zenith angle and energy bins. The expected 
SSM flux of a particular energy bin is calculated from 
the total and hep flux of BP2000|§ and the neutrino 
spectrum from Ortiz et a/. Jl^. This neutrino spectrum 
is based on an improved measurement of the /3-delayed a 
spectrum of the *B decay with a small and well-controlled 
systematic uncertainty. Earlier reports ||l^ used the neu- 
trino spectrum by Bahcall et al. p^] . 

This zenith angle spectrum is analyzed in a two- 
neutrino oscillation scenario, which can be described with 
a mixing angle and a mass difference Am^. We con- 
sider two cases: (i) Ve <-> V)j.,t and (ii) Ue ^ I'stGiiic- For 
each set of neutrino oscillation parameters (sin^ 26 and 
Am?) the expected number of solar neutrinos and its 
zenith angle spectrum are calculated. First, the proba- 
bility Pi (P2) of a solar neutrino to be in the mass eigen- 
state vi (1^2) on the surface of the sun is obtained from a 
numerical calculation which propagates a neutrino wave 
function from the production point in the core to the sur- 
face. This calculation uses models for the distribution of 
the neutrino production point in the sun , the electron 
density in the sun and the neutrino spectrum [ p^ . 
Above Am^ — 1.8 • 10~^ eV^ the propagation of the two 
mass eigenstates from the sun to the earth and inside 
the earth can be assumed to be incoherent. The survival 
probability at the detector is given by 

P{l^e ^ '^e)sK = PlPle + P2P2e, (1) 

where Pie (P2e) is the probability to be v^, at the detec- 
tor if the neutrino arrives at the earth as vi {V2) taking 
into account matter effects inside the earth pi. The elec- 
tron density model for the earth (PREM |14|) assumes 
a charge-to- mass ratio (Z/A) of 0.468 for the core and 
0.497 for the mantle ||l||. Below Am^ = i.g . IQ-^ eV^ 
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eV^). The SMA and just-so solutions predict spectral 
distortion, while the LMA and LOW solutions predict a 
zenith angle variation and can therefore be constrained 
by the zenith angle spectrum data of SK. The hatched ar- 
eas of Fig. ||are the combined analysis allowed regions at 
95% C.L. assuming i^e ^ '^sterile oscillation. In this case, 
LMA and LOW solutions do not occur, since there is not 
enough neutral current contribution to neutrino-electron 
scattering to accommodate the difference between the SK 
and the Homestake flux results. 

Using the zenith angle spectrum (see Fig. |l|) the prob- 
ability of two neutrino oscillation scenarios was tested 
with a method. For each energy bin i, we form a 
zenith angle flux difference vector A^. Its seven zenith 
components A^^^ are 



A. , = 



a X f {Ei, (5<,orr) X 



where z is the zenith angle bin, 0"°"° is the observed flux 
of each energy and zenith angle bin, and ^f^" and (j)""^ 
are the expected event rates in that bin without and with 
neutrino oscillation. The spectral distortion / due to the 
correlated systematic error (see table |) of ipi^z is scaled 
by the parameter S^arr- The definition of the is thus 



FIG. 1: Spectrum between 5 and 20 MeV for various zenith 
angles. Data points with open circles combine all zenith angle 
bins. The error bar is the statistical error. See table | for 
systematic errors and the definition of the binnig. 



matter effects inside the earth are unimportant and the 
propagation from the sun to the earth is assumed to be 
coherent. The survival probability is then 

P{i^e ^ '^e)sK = Pi cos^ 6* -f P2 sin^ 6*+ 

2 X v/PiP2 cos 61 sin 6>cos( ^"'' ^ )i2) 

where L is the distance from the sun to the earth ranging 
from perihelion (winter) to aphelion (summer) [ pO| . 

Four experiments have quoted measurements of the 
solar neutrino flux: Homestake (2.56 ± 0.16 ± 0.16 
SNU), SAGE § (75.4+^j SNU), GALLEX/GNO § 
(74.1+^;^ SNU), and SK g. Taking an average of SAGE 



and GALLEX, a combined analysis results in several al- 
lowed regions. The combined analysis is performed by 
the method given in |^ considering updated theoretical 
correlated and uncorrelated uncertainties ("^^Cl cross sec- 
tion |l^, ^^Ga cross section ||l^, neutrino-electron scat- 
tering cross section and diffusion jl^). The hatched 
areas of Fig. || show the allowed regions at 95% C.L.(x^ < 



5.99) assuming Vf, 



oscillation. There are 



four allowed regions called "small mixing angle solution" 
(SMA, Ato^ « 10-5 eV^, sin^ 29 « IQ-^ . . . 10"^), "large 
mixing angle solution" (LMA, Ato^ ~ iQ-* . . . iQ-5 gy^, 
sm^2e > 0.5), "low solution" (LOW, Am^ « lO"'^ eV^, 
sin^26l « 0.9) pi, and "just-so solution" (Am^ < lO'^ 



i=l 



(3) 



Each energy bin i has also a separate 7x7 error matrix 
Vi describing the energy-uncorrelated uncertainty. Vi is 
the sum of the statistical error matrix and the energy- 
uncorrelated systematic error matrix (see table ^, the 
latter of which is constructed assuming full correlation 
in zenith angle. The flux normalization factor a is un- 
constrained to make the independent of the total so- 
lar neutrino flux. The correlation parameter S^an- is con- 
strained within cr^o„ . The size and shape of the correlated 
error are calculated as in [0. The hep contribution to 
the neutrino flux is not constrained. 

The values are calculated in the parameter space, 
(lO"'^ < sin^ 20 < 1, 10"" eV^ < Am^ < lO'^ eV^). In 
the case of active neutrinos, the minimum x^ value is 36.1 
with 40 degrees of freedom at (sin'^ 26 = 1, Am^ — 6.53 • 
10^1" eV^). The best-fit flux normalization is a = 0.788, 
the correlation parameter is 6^„,. = ~0.06a^„, and the 
hep flux is 0. The shaded areas in Fig. || are excluded at 
95% C.L. from this flux independent analysis. Most of 
the SMA and just-so solutions are disfavored with this 
C.L.. In the case of sterile neutrinos, the minimum x^ 
value is 35.7 at (sin^ 20 = 1, Am^ = 6.57 ■ lO^^" eV^). 
All possible solutions are disfavored at 95% C.L. in this 
case. The best-fit fiux normalization is a = 0.917, the 
correlation parameter is S^^, — 0.06(7^a„ and the hep fiux 
is 0. The shaded areas in Fig. ^ show the excluded regions 
(95% C.L.). 

Using the theoretical uncertainty of the fiux 
ffiux =lo 16^^^' analysis combining flux and zenith 
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FIG. 2; Exclusion area for two-flavor osciUation ^ fp/i^r 
from zenith angle spectrum analysis at 95% confidence level. 
Overlaid are the allowed areas (95% C.L.) using the zenith 
angle spectrum and the SSM flux prediction (dotted lines). 
The small overlap of allowed area and excluded area is due 
to the fairly good agreement of the flux for these parameters. 
The hatched areas are allowed at 95% in a combined flt to 
the fluxes measured at GALLEX§, SAGE|||, Homestake^] 
and Super-Kamiokandep]. 



angle spectrum has also been performed. In the active 
neutrino case, the minimum value is 37.8 with 41 d.o.f. 
at the same position as the unconstrained case. The 
flux normalization changes to a = 0.789 and the cor- 
relation parameter to = — 0.02tTcorr The minimum 
point is within the just-so solution, but some LMA 
are similar to the minimum. For example, ^ = 39.1 
at (sin^ 26* = 0.87, Am^ = 7 • IQ-^ eV^) with a hev flux 
of 2.9xBP2000. The dotted lines in Fig. |show the con- 
tours of the 95% C.L. allowed regions. In the sterile neu- 
trino case, the minimum value is 35.9 with 41 d.o.f at 
the same position as the unconstrained case. Flux nor- 
malization, correlation parameter and /lep flux are un- 
changed. The inside of dotted lines in Fig. || is allowed 
at 95% C.L.. Since the allowed area from the combined 
flux analysis does not overlap these regions, oscillations 
into only sterile neutrinos are disfavored at this confi- 
dence level. 

Figures |^ and ^ are based on the analysis of the 



FIG. 3: Exclusion area for two-flavor oscillation ^ J^storile 
from zenith angle spectrum analysis at 95% confldence level. 
Overlaid are the allowed areas (95% C.L.) using the zenith 
angle spectrum and the SSM flux prediction (dotted lines). 
The hatched areas are allowed at 95% in a combined flt to the 
fluxes measured at GALLEX|| and SAGEg, Homestake^] 
and Super-Kamiokandep]. 



zenith angle spectrum. We have also performed an oscil- 
lation search using the "day/night spectrum", which, in 
contrast to the zenith angle spectrum, divides the data 
into two zenith angle bins (day and night bin). Each of 
these bins is then divided into 19 energy bins The 
is defined as follows: 

19 

D,JV 1=1 



The notation is analogous to that used in the defini- 
tion of the zenith angle spectrum analysis, ai is the sum 
of statistical and uncorrelated errors added in quadra- 
ture. 

The minimum y^ value is 28.2 with 34 degrees of free- 
dom at sin^ 29 = 0.4 and Am^ = 1.38 • 10"!° eV^. The 
best-fit flux normalization is a = 0.488 and the correla- 
tion parameter is = — 0.2(Tco„. Figure shows the 



5 



-3 

cv, 10 
> 

(L> 




-8 

10 r 



-9 

10 r 




-4 -3 -2 -1 

10 10 10 10 , 1 

sin^(20) 



FIG. 4: Exclusion area for two-flavor oscillation v^. <-> v^jvr 
from a day/night spectrum analysis. The hatched area is 
allowed at 95%C.L. by the combined flux analysis same as in 
Fig. I 

95% excluded regions using the shape of this day/night 
spectrum. The excluded area is similar to that obtained 
in the zenith angle spectrum, but more restrictive in 



the SMA region. The differences at the LMA and near 
the LOW solution are due to the zenith angle variations 
within the night bin. The lower left corner of the SMA 
predicts a slight depression of the core flux resulting in 
a day flux prediction that is larger than the night flux. 
SK measures a l.Scr excess of the night flux over the day 
fluxQ, but the flux in the core bin is below the day flux. 
This leads to a slightly better fit of these parameters to 
the zenith angle spectrum than to the day/night spec- 
trum. The lower left corner of the SMA 95% C.L. region 
is excluded at 93% C.L by the zenith angle spectrum and 
at 97% C.L. by the day/night spectrum analysis. Other 
differences are due to the use of different binnings. 

In summary, Super-Kamiokande precisely measured 
the energy dependence and zenith angle dependence of 
the solar ®B neutrino flux. The data do not show a sig- 
nificant distortion of the spectrum or zenith angle vari- 
ation. This places strong constrainins on neutrino os- 
cillation solutions to the solar neutrino problem inde- 
pendently of the flux expectation. If oscillations into 
active neutrinos are assumed, just-so and the SMA so- 
lutions are disfavored at 93% (zenith angle spectrum) to 
97% C.L. (day/night spectrum) and the LMA solutions 
are preferred. In conjunction with the SK flux mea- 
surement, two allowed areas at large mixing remain. All 
possible oscillation solutions into only sterile neutrinos 
are disfavored at 95% confidence level. 
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